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ABSTRACT 


The lunar electrical conductivity profile has been 
obtained from an iterative fit of theoretical electromagnetic 
transfer functions to empirical ones based on joint power 
spectral density analyses of data from the Apollo 12 Lunar 
Surface Magnetometer and the Ames Explorer 35 magnetometer. 

Seven selected two hour swaths and seven selected one hour 
swaths of data were used. The spectrum analyzed ranged in 
frequency from 0.00083 to 0.04 Hz. The amplification of the 
interplanetary magnetic field at the lunar surface showed a 
distinct increase from ~ 1 at 0.001 Hz to ~ 3 at 0.005 Hz. 

At higher frequencies the slope of the amplification vs. 
frequency curve decreases until the amplification levels off 
at about 4 at 0.025 Hz. An ampl ification curve with such a 
distinct bend is best fit by a lunar conductivity model with 
a sharp maximum around 1500 km radius, so that high frequency 
variations in the magnetic field are compressed into the outer 
shell while low frequency variations can penetrate to the 
interior. For a monotonic temperature profile, the electrical 
conductivity model requi res a change in composition around 
1450 km radius. A thermal and compositional model which appears 
to fit the data is a basalt-like outer layer, an olivine-like 
core and a temperature of 450°C at the conductivity peak, 
increasing to 800°C in the deep core. 
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INTRODUCTION 

This paper extends our preliminary report (SONETT et al . , 
1971a) of the bulk electrical conductivity profile of the Moon 
determined by analysis of data from the Apollo 12 Lunar Surface 
Magnetometer (LSM) and the Ames Explorer 35 lunar orbiter 
magnetometer. The account given here is still exploratory but 
more details are supplied on the mechanism of induction, the 
analysis, and the implications of our findings. 

It has been recognized in recent years that electro- 
magnetic excitation of the Moon by signals arising in the 
solar wind produce electrical currents in the deep lunar 
interior. The magnetic fields associated with these currents 
are detectable at the surface of the Moon (SONETT et al., 1971a, 
b,c; DYAL et al., 1970 ) and, in the case of a strong 
inductive response may even modify the plasma flow near the 
Moon. The establishment of Explorer 35 in lunar or bit permitted 
study of the plasma environment close to the Moon. However, 
the detection of electromagnetic events in the Moon appears to 
be feasible only from either the surface or a very low altitude 
orbiter. Using the LSM, lunar electromagnetic induction can be 
assured. We believe the corresponding theory is sufficiently 
well understood that lunar magnetometry can be used to determine 
the electrical conductivity of the interior and thence to make 
inferences of temperature and composition. Although a detailed 
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analysis of spatial resolution has so far not been carried 
out, our preliminary assessment is that volume elements of the 
order of 3% are resolvable in the region of the steep con- 
ductivity rise with depth. 

A fundamental difference of the Moon from the earth arises 
from the solar wind. When the angle between the Moon-sun line 
and the local vertical at the LSM is appreciably less than 90°, 
the dynamic pressure of the solar wind is known to compress 
the induced magnetic field into the less conducting subsurface 
layers of the Moon (SONETT et al., 1971a), appreciably amplify- 
ing the induction over that in a pure vacuum or a nonconducting 
atmosphere. This distinction together with the presence of 
the diamagnetic cavity results in a day-night induction asymmetry. 

MECHANISM OF THE INDUCTION 

Electromagnetic induction in the Moon is dependent upon 
a variety of discontinuities and waves in the solar wind. The 
forcing function which drives the induction is composed of the 
steady interplanetary magnetic field upon which is superimposed 
a hydromagnet i c radiation continuum due to ( 1 ) plasma waves 
arising presumably in the solar atmosphere and subsequently 
convected and propagated outwards, (2) waves due to local 
instabilities in the solar wind, and (3) discrete large amplitude 
events such as collision-free shock waves, tangential dis- 
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continuities, and Alfv6n waves (cf . COLBURN and SONETT, 1966"). 

The electromagnetic interaction includes toroidal and 
poloidal magnetic fields. The principle excitations for these 
are a) the interplanetary electric field (in a reference 
frame comoving with the Moon) given by _E m = v x _B . , where v 
is the velocity of the Moon relative to the solar wind and _B 

the instantaneous interplanetary magnetic field, and b) the 

« 

time rate of change of the interplanetary field, J3 . The two 
modes correspond respectively to transverse magnetic (TM) and 
transverse electric (TE) excitation. Both modes display strong 
frequency (f) dependence. The TE transfer function is zero at 
f = 0 and increases with increasing f . The TE mode currents 
which close wholly in the lunar interior, tend to be concentrated 
where variations in the magnetic field are damped substantially. 
With increasing frequency, the poloidal fields become compressed 
into shells of decreasing thickness, thus increasing the amplifi- 
cation of the magnetic field as observed at the Moon's surface. 
Examination of the interior of the Moon to depths around 800 km 
is facilitated by this mode, but for greater depths is limited 
by the decreased response. 

The TM mode attains peak value for f = 0 ; it remains 
approximately constant with increasing f until a combination 
of core and lithospheric conductivities forces the currents to 

* Although a large body of experimental work has subsequently 
been published, this paper provides the essential background 
for large events. 
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pass wholly through the crust whereupon it decreases with 

further increase in f beyond about 0.03 Hz. This mode is 

responsible for steady state bow wave phenomena. Attempts to 

detect a bow wave in the solar wind upstream from the Moon 

have not met with success. The TM current system, which is 

required to pass through the lunar crust and close in the solar 

q 

wind, is consequently limited to less than 10 amperes. How- 
ever, as shown later, the total absence of the TM mode is not 
assured and evidence for this mode is present. 

Theoretical treatments of the lunar interaction with the 
solar wind generally include the effect of the solar wind 
dynamic pressure in confining the induced field lines (SONETT 
and COLBURN, 1968; JOHNSON and MIDGLEY, 1968; BLANK and SILL, 
1969; SCHUBERT and SCHWARTZ, 1969). This effect is provided 
in the model by a field confining surface current layer in 
the solar wind just ahead of the lunar surface. Preliminary 
examination of the lunar response using the LSM data shows a 
strong amplification of incident tangential discontinuities 
whose free stream properties are monitored by Explorer 35. 

The amplification occurs only for the vector components 
tangential to the surface; the normal component tends to follow 
the interplanetary value. Thus the existence of a confining 
current layer appears verified for the sunward side of the Moon. 
The very strong amplification implies that the lines of force 
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are confined within the Moon to a crustal layer having an 
electrical conductivity substantially less than that of the 
deeper 1 ayers( SONETT et al . , 1971b). 

THEORY OF THE INDUCTION 

A complete treatment of the electromagnetic interaction 
would requi re appropriate matching of the interior fields to 
those in the plasma surrounding the Moon. The exterior flow 
field shows no marked perturbations aside from the diamagnetic 
cavity. An analysis based upon an inhomogeneous Moon immersed 
within a perfectly conducting space is used for the sunward 
hemisphere. This simplification is justified since cavity 
current effects are much smaller on the sunward hemisphere. 

We consider the interaction of a solar wind forcing 
magnetic field oscillation 

H = j] Hg exp|2TT i (£ - ft)} , ( 1 ) 

with a radially inhomogeneous Moon. The cartesian coordinate 

A A A 

system (?»ri,C) with unit vectors is fixed relative to 

the Moon and has its origin at the Moon's center. The Moon 
moves with speed v in the negative £ direction. The 
quantities Hq and X = v/f are the amplitude and wavelength, 
respectively, of the magnetic field oscillation. The magnetic 
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field forcing function alone drives the TE lunar response. A 
derivation for the TM mode response is not included here be- 
cause it must be relatively unimportant (see the discussion 
in the following section). In the lunar interior the solution 
of Maxwell's equations for the TE mode can be represented by 
the potential Q which satisfies 

v 2 Q + k 2 Q = 0 , ( 2) 


where 


2 2 

k = uo pe + iapuu , uu = 2rrf , (3) 

a is the electrical conductivity, and p and c are 
respectively the magnetic permeability and permittivity of 
the Moon. Although some possibility exists that local effects 
are significant in increasing the value of p over the free 
space value, there is presently no evidence for this. DYAL 
and PARKIN (1971) give (1.03 ± 0.03)Pq for the global perme- 
ability. In the following discussion the values assumed 
for the global permeability and permittivity of the Moon 
are the free space values of these quantities. Possible 
departures of the permittivity from the free space value 
are insignificant in the subsequent application of this 
theory. The equations for determining H from the potential 
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Q are presented in the Appendix. 

The boundary condition applied to the TE mode is the 
continuity of the normal component of the magnetic field at 
the lunar surface. This condition results from the confine- 
ment of the induced field, treated as a current sheet at the 
Moon-plasma interface. 


The potential 0 is given by 


0 = pvHq y sincp ^ 3^G^( r)pj(cos0 ) 

1=1 


(4) 


where (r,9,cp) are spherical polar coordinates with C the 




( 21 + 1 ) 


polar axis, a is the lunar radius, 3^ is £+'] j > 
p] are associated Legendre polynomials. The functions G„ 
are solutions of 


and 


d 2 G 


t. . J,.2 


+ k 


dr‘ 


{' 


g 4 = o 


(5) 


with the boundary conditions 


<=* < r = a > ■ h t-r 2 ) 


(6) 


where are the spherical Bessel functions. The equation 

for the potential and the form of the boundary conditions 
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follow, in part, from the character of the spherical harmonic 
expansion of the solar wind forcing field, given in the 
Appendix. LAHIRI and PRICE (1939) were the first to obtain 
the radial differential equation for the TE mode in connection 
with their investigation of the geomagnetic induction problem. 


For the interpretation of LSM data we introduce the modal 
transfer function for the tangential components of the TE 
mode magnetic field at the lunar surface 



( see Append ix) 


INSTRUMENTATION AND DATA 

General properties of both the LSM and the Ames magnetom- 
eter on Explorer 35 are given elsewhere ( DYAL et al . , 1970; 
MIHALOV et al . , 1968). Provision has been made in the Explorer 
magnetometer for the suppression of spin tone modulation of the 
data spectrum by utilizing a pair of synchronous demodulators 
which operate in quadrature upon the two spin tone modulated 
signals. This suppression is done on the spacecraft and can 
be shown to lead to a time series from which spin modulation is 
ideally eliminated. The signals returned to Earth are formally 
identical to those obtained from a nonspinning inertial ly stable 
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spacecraft (SONETT, 1966). 

A sample and hold system samples the three components of 
the vector at a uniform rate. In order to avoid the risk of 
unacceptable alias the magnetometer outputs are passed through 
a low pass filter prior to sampling, according to the Nyquist 
criterion. The effect of filter rolloff is significant for 
frequencies above 0.025 Hz, which is the 3 db point. The power 
spectral densities have been corrected for this effect at all 
f requenci es . 

Spectra have been obtained both for cases where the Moon 
is in the magnetosheath and also in the free stream solar wind. 
The forcing field defined by Explorer 35 measurements is trans- 
formed into the local LSM coordinate system (x is along the 
normal outwards from the surface, y is easterly and z is 
northerly at the site of ALSEP and the LSM). The LSM data, 
which define the response field, are edited to eliminate 
spurious data and to insure time continuity. Time gaps and noise 
in the LSM data are aperiodic and are attributed to telemetry 
transmission and data processing. For the cases studied the 
amount of data missing or deleted because of gaps and spurious 
noise is small (< 5%). The data are then numerically filtered 
and decimated in order to analyze segments of the record with 
approximately the same sampling interval and upper frequency 
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limit as that from Explorer 35, the latter having a Nyquist 
limit approximately one order less than LSM. The Explorer data 
are similarly edited to insure time continuity and to remove 
obvious spurious signals. For the Explorer data the time series 
are continuous for the intervals studied. Spurious signals at 
random intervals, are removed and again constitute a small 
amount of the time interval (< 5%). 

The power spectral densities for the forcing and response 
fields are computed using standard techniques for determining 
the autocorrelation function, smoothing, and taking the Fourier 
transform (BENDAT and PIERSOL, 1966; JENKINS and WATTS, 1968). 
The spectra are analyzed using 20 degrees of freedom, and the 
autocorrelation function is smoothed with a Parzen weighting 
function (JENKINS and WATTS, 1968). This gives an approximate 
error of 25 per cent in each power spectral density estimate. 

We have applied these techniques to seven selected two hour 
swaths and seven selected one hour swaths of data. Longer 
swaths, which would have reduced the horizontal error bars in 
Fig. 2, were not used because of time gaps, on the order of 
10-30 minutes, in either the LSM or Explorer data. These are 
caused by telemetry shadowing, calibration interruptions, etc. 

A representative set of spectra at both Explorer and LSM 
for the y and z magnetic field components (tangent to the 
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surface) is shown in Fig. 1, The forcing spectra are seen to 
display the expected f dependence character! sti c of inter- 
planetary fluctuations, while the LSM shows an f"^ dependence. 
Thus significant frequency dependent power amplification is 
apparent in the record. 

The empirical transfer function A j ( f ) is defined by 

h 2i (f) + h n (f) = A, (f) h n (f) , (7) 

where (f) and h-j j (f) are the Fourier transformed time 

series of the magnetic field induced in the Moon and measured 
on the lunar surface, and the free stream interplanetary 
magnetic field, respectively, and the subscript i is x, y 
or z. The LSM measures the sum (f) + h^ j (f) while the 

free stream magnetic field, h^ . (f) , is measured by the Ames 
magnetometer on the lunar orbiting satellite, Explorer 35. 

Fig, 2 shows the average ratio of power spectral densities 
Aj(f) . The mean values are the arithmetic averages using the 
14 swaths of data. The error bars are the one standard devia- 
tion limits of the mean calculated from the transfer functions 
for the 14 time swaths. The bars do not include the error 
estimates in the calculation of the individual spectra. They 
are, however, consistent with the propagation of errors in the 
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individual spectra. The values of A lie near unity at all 

y\ 

frequencies indicating that the interior conducting region is 
relatively far from the magnetometer compared to the distance 
to the plasma current sheath on the sunward side of the Moon. 
Characteristically the amplification rises from near unity for 
both Ay and A z at 0.001 Hz until values near 4 are attained 
above about 0.02 Hz. Values of A y are less than those of 
A z over the entire frequency range. This departure increases 
fractionally as frequency is decreased. It is plausible that 
the lunar TM response is frequency dependent at frequencies 
above about 0.01 Hz. If so the response for the TM mode should 
decrease with frequency in accordance with the observation. 

This is the strongest evidence for contributions from both the 
TE and TM modes. An alternate possibility that the basic 
response is purely TE and that the Moon is asymmetrically ex- 
cited seems less convincing but cannot be ruled out at present. 
We have not yet found a completely satisfactory explanation for 
why the northward component of the excitation is systematically 
higher than the eastward component. The important facts are 
that they both have the same slope and change in slope around 
0.01 Hz. We conclude that the lunar transfer function is 
dominated by the TE mode with possible contributions, especially 
at the lower frequencies, from the TM mode. 

A significant scatter appears in the amplifications Aj . 
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This scatter is substantially reduced by taking averages over 
many Aj , but is exact only if the errors have zero mean, 
which is approximately the case as shown in Fig. 3. A combina- 
tion of all the data points "A(f) = (A^(f) + A^(f))j- 2 for 

all frequencies up to 0.035 Hz has been made in order to 
investigate the distribution of the ratios of the power 
spectral densities. The distribution shown in Fig. 3 is a 
histogram of the number of cases vs. the difference of the 
measured value “A from the mean at each frequency, normalized 
by the standard deviation at that frequency and weighted by 
the expected error ( ^ ( f q ) ( f ) ) • The figure shows that the 
distribution is somewhat more sharply peaked than a normal one 
and also has a small bias. 

The sources of the variations in the A. may originate in 
any of several possible ways: 1) the presence of contamination 

especially in the Explorer data from spurious tones in the 
spacecraft; 2) the contribution to induction from higher order 
modes; 3) fluctuations in cavity currents in turn due to 
changing diamagnetism in the solar wind; 4) the possible contri 
bution of TM modes to the predominantly TE spectra; 5) the poss 
bility that some field line leakage takes place where the 
component of solar wind plasma pressure normal to the surface 
falls below the pressure of the induced magnetic field. These 
factors may also contribute to the droop of the average Aj at 
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the higher frequencies. It should be noted that there is no 
correlation between the Aj and sun angle, implying that field 
line confinement is effective for the data used here (sun angles 
1 ess than 65° ) . 

LUNAR ELECTRICAL CONDUCTIVITY PROFILE 

The theory of lunar induction, discussed earlier, has 
been used to derive a profile of electrical conductivity from 
the empirical transfer function. The theoretical amplification 
for the lowest TE mode, computed by numerically integrating 
Eq. (5) with a radially dependent conductivity, is matched to 
the empirical transfer function using a Newton- Raphson iterative 
scheme. This scheme readjusts the entire conductivity profile 
to yield amplifications which best fit the data in a least- 
squares sense over the entire frequency range. This has been 
carried out using frequency values of 0.83, 1.7, 5, 12, 17, 22, 
25 and 35 millihz. The conductivity profile is characterized 
by its values at the radial locations r = 800, 1200, 1400, 

1450, 1490, 1510, 1550 and 1740 km. For r < 800 km the con- 
ductivity a is set to the value at r = 800 km ; elsewhere 
a linear interpolation of log a is used. 

The computer calculation is started with a continuous 
conductivity profile defined by 8 parameters as described 
above. Numerical integrations are carried out to obtain 
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values of amplification at the eight frequencies, A comparison 
of these amplifications with the data provides the Newton- 
Raphson scheme with the input required to adjust the 8 con- 
ductivity parameters to yield a profile whose corresponding 
theoretical amplification curve is a better 1 east- squa res fit 
to the empirical transfer function. The iteration is continued 
until an adequate fit is obtained. The result of such a fit 
to the empirical amplification curve = {^-(Ay + A^)} 2 is 
shown in Fig. 4. The introduction of 7? permits the use of 
both Ay and A z data to improve the statistics. The dif- 
ferences between the theoretical amplification and the empiri- 
cal 7\ are, in our assessment, partially attributable to the 
various complications of the excitation process enumerated at 
the end of the previous section (not accounted for by the 
theory), as well as considerations of computer running time 
which limits the number of frequencies and conductivity para- 
meters (spatial resolution) employed. 

A conductivity profile is associated with the theoretical 
amplification curve of Fig. 4. In addition we have obtained 
conductivity profiles by fitting not only 7C but also Ay , 

A z and their one standard deviation limits. The conductivity 
profiles derived from the iterative least squares inversion 
for Ay , A z and J\ are shown in Fig. 5. The prominent spike 
in each conductivity profile is an invariant characteristic 
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of the inversions. It is centered at about r = 1500 km where 
the conductivity is nearly 10 mhos/m. The inner minimum lies 
at about r = 1400 km and the conductivity appears to rise at 
greater depth. Gross bounds on the conductivity profile are 
seen in the insert; these are determined from the one standard 
deviation limits of the various A's, but do not themselves 
represent one standard deviation limits on the conductivity 
prof i 1 e . 

The computer calculations for Figs. 4 and 5 started with 
a constant conductivity of 1 0 ^ mhos/m. However, a number of 
computations have been carried out using different values for 
uniform starting conductivities ( e . g . 10 J and 10 3 mhos/m) and 
different radial locations. In every case tested, initial 
convergence was rapid and the final conductivity profile invari- 
ably displayed the prominent spike near r = 1500 km. For the 
profiles reported here, several values of r were chosen in 
the neighborhood of r = 1500 km to better define the con- 
ductivity spike. 

Whereas the large spike in conductivity is a persistent 
feature of the inversions, the character of the conductivity 
profile at greater depth, where the conductivity appears to 
rise, is not so certain and our results for the conductivity 
at these depths must remain tentative. The surface amplifica- 
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tion Is a rather insensitive function of core conductivity 
because of the relatively small core volume and the distance 
to the surface. 

Figure 6 shows the manner in which the amplitude and 
phase of the relative field strength vary with depth in the 
Moon, indicating different degrees of field compression at 
various frequencies. The conductivity profile used in this 
calculation was the one associated with the theoretical fit 
of (see Fig. 5). In the range r = 1500 to 1 7^+0 km , the 
field is essentially constant and phase shift is small, con- 
sistent with the low conductivity (large skin depth) in this 
region. The large rise in conductivity at r = 1500 km 
produces a penetration barrier (small skin depth) for the 
higher frequencies. Below this barrier the high frequency 
variations in the field are low in amplitude and have large 
phase shifts. Over the high frequency range the barrier 
essentially confines field lines to the outer region, produc- 
ing a nearly constant amplification with frequency, consistent 
with the data. For the lower frequencies the barrier is not 
effective and the relatively low conductivity in the range 
r < 1500 km allows the field to penetrate to the deep interior, 
yielding a smaller amplification. Where phase shifts are greater 
than 90 deg the field lines connected with. the driving field 
have been excluded; the resulting compression leads to the 
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amp I i f f cat i on . 

COMPARISON WITH OTHER CONDUCTIVITY PROFILES 

A physical understanding of the appearance of the large 
conductivity spike can be obtained as follows. In a two 
layer model with an infinitely conducting core and a non- 
conducting shell the amplification of the tangential magnetic 
field components is given by 

, ,3 (core volume) 

1 + 1 (shell volume) 

This simple result requires that the high frequency skin depth 
barrier be near r = 1500 km for an amplification « 4 . At 
frequencies above 0.02 Hz essentially no wave, penet rat i on takes 
place through this barrier. Thus the volume available in the 
core for field lines is insignificant compared to that in the 
nonconducting shell. This means that the amplification be- 
comes independent of frequency, i.e. dA/df -* 0 at high 
frequency. Concurrent with the high frequency limitation is 
the requirement that A drop to near unity at the lower 
frequencies. If it were assumed that the electrical con- 
ductivity were monotonic, so that the interior conductivity 
were uniformly high, then the low frequency amplification 
would be in excess of the observed values. 
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The inadequacies of monotonic conductivity profiles are 
shown in Fig. 1. The experimental r.m.s. amplification data 
A (circled points with associated error bars) can be compared 
with the amplification curve for the A conductivity profile 
of Fig. 5 (labelled "best fit") and amplification curves of 
several monotonic conductivity profiles. The "2 layer" ampli- 
fication curve was calculated for a model with a core of radius 
of 1560 km and a constant conductivity of 7.6x10”^ mhos/m, and 
a shell of zero conductivity. This value of the core con- 
ductivity provides a best fit to the experimental A for the 
given core radius. Other best fit two layer models with 
different core radii have been investigated; the one used in 
Fig. 7 yields an amplification curve which matches the data 
most closely. The amplification curve labelled DYAL and PARKIN 

(1971 ) was computed from a 3 layer conductivity model proposed 

- 2 

by these authors. The parameters of their model are a - 10 
mhos/m for 0 < r < 1 044 km, a = 1.7xl0 -if mhos/m for 1044 km < 
r < 1653 km and a = 0 in the outer shell. In other 3 layer 
conductivity models consistent with the data analysis of DYAL 
and PARKIN (1971), the outer boundary of the region of inter- 
mediate conductivity can be located anywhere between r = 1653 
and 1740 km. We have computed the amplification curves for a 
number of these additional 3 layer models; none provide a 
better fit to the observational curve. Also shown in Fig. 7 
are amplification curves from conductivity profiles proposed 
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by SILL (1971) (model 3 in that paper) and NESS (1969). 

Other conductivity models investigated by SILL (1971) fit the 
experimental data no better than his model 3« The conductivity 
model of NESS (1969) consists of a core of radius 1426 km 
with o = 8x10 ^ mhos/m, and a nonconducting shell. Amplifi- 
cation curves from conductivity models of WARD (1969) were 
also considered but these were rejected upon comparison with 
the observational data. 

The "2 layer" amplification curve shown in Fig. 7 is, 
compared to the experimental data, high at f < 0.0035 Hz and 
f > 0.03 Hz, and low in the range 0.01 Hz < f < 0.02 Hz. The 
data show a flat response at frequencies above 0.02 Hz, whereas 
the slope of the "2 layer" curve is high at these frequencies. 
Furthermore as the frequency increases the curvature of the 
data changes from positive to negative at about 0.003 Hz. The 
"2 layer" ampl i f i cat i on curve is everywhere concave up. Th i s 
"s" shaped character of the experimental data is faithfully 
reproduced only by our "best fit" conductivity profile . Further 
work is in progress to improve the high frequency slope of our 
"best fit" model. This will probably lead both to a higher 
value of the maximum conductivity and a steeper slope of the 
conductivity profile at the outer edge, effects which will 
tend to sharpen the conductivity spike in order to allow 
penetration of low frequency magnetic field fluctuations. 
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COMPOSITIONAL AND THERMAL MODEL 

The electrical conductivity profile cannot be explained 
by a uniform material and a plausible thermal profile. In 
the region from the surface to r = 1500 km where the con- 
ductivity attains its maximum value, the rise of conductivity 
with depth is a reasonable consequence of the accompanying 
increase of temperature in a material of uniform composition. 
Below r = 1400 km the apparent rise in conductivity is again 
explainable by an increase in temperature. On the other hand 
the precipitous decrease of electrical conductivity by 2 to 3 
orders between r = 1500 and 1400 km cannot be explained as 
due to temperature. Either a compositional change, phase 
change, or a combination of the two is required. Thus a 
reasonable model for stratification of the Moon, limited by 
the present poor spatial resolution of the analysis indicates 
a core out to r « 1400 km overlain by a mantle of higher 
conductivity material, plus possibly a transition layer at 
r = 1 400 - 1500 km. 

In order to infer a temperature profile from the con- 
ductivity profile it is necessary to use conduct i vi ty- tempera- 
ture functions of known rock materials. The Apollo basalts 
( SCHWERER et al., 1970; NAGATA et al . , 1970) are representative 
of the most conducting rocky matter known. On the other hand 
olivine (dunite) or ol i vi ne- per i dot i te (ENGLAND et al., 1968) 
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is representative of poorly conducting geological material. 

The peak conductivity value found in the lunar mantle corre- 
sponds to a temperature of about 450°C for lunar basalt or 
950°C for olivine. Since the temperature just under the lunar 
surface is -30°C, the corresponding thermal gradients in the 
outer mantle are about 2°K/km and 4°K/km. 

The precipitous decrease in electrical conductivity of 
2 3 

10 - 10 as depth increases from 250 km to 350 km takes place 

in a distance where the temperature rises at least 100 degrees. 
While the absolute conductivities of the materials are unsure 
because of impurities, etc., it seems clear that the change in 
conductivity around 300 km depth requires a change in composition. 
Thus one model fitting the data has an olivine-like core, 
possibly a lower mantle composed of a material transitional 
between the substrate and the upper mantle, and an upper 
mantle having the conductivity of Apollo basalt (see Fig. 8). 

Since the present frequency resolution is insufficient to 
estimate the conductivity close to the surface, these results 
do not preclude the presence of an anorthositic crust (WOOD 
et al„, 1970). The volume of the mantle comprises about one 
half the total volume of the Moon. To f i 1 1 it completely with 
a basaltic extract is in conflict with geochemical constraints 
which imply at most a 10-15% contribution by volume. The upper 
mantle may have to be composed only partly of a basalt- like 
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substance to account for the apparent dependence of the upper 
mantle conductivity upon temperature. The actual basalt-like 
thickness cannot be assayed yet because of limited spatial 
resolution. An alternative model which avoids the difficulty 
of finding a source of so much basalt has a layer of iron con- 
centration at r = 1400 to 1500 km, such as could occur from 
a differentiation of the outer 340 km (WOOD et al., 1970; cf . 
UREY et al . , 1971 ). 

For our compositional model the temperature at r = 1500 km 
Is 450°C and the average near surface thermal gradient is 
2°K/km. The temperature at r = 800 km, computed from the con- 
ductivity profiles shown in Fig, 5 together with an olivine 
conductivity function is in the range 750-800°C, decreasing 
to approximately 550°C at 1400 km where the local minimum in 
electrical conductivity is observed, A calculation based upon 
more conducting matter would depress the computed temperatures. 
R* REYNOLDS (private communication) has shown that after 4,5 
billion years a nonconvect i ng Moon of 25% ordinary chondritic 
radioactive composition and an initial uniform temperature of 
0°C attains a near surface thermal gradient of approximately 
2°K/km and a uniform core temperature of about 800°C, Hence 
a thermal model consistent with the conductivity requires 
appreciable depletion of the heat sources in the Moon relative 
to ordinary chondrites. 
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An approximate value for the lunar heat flux can be cal- 
culated from the near surface gradient determined above. For 

silicate material the thermal conductivity is 0.008 cal/cm sec 
deg though it is known that it will decrease with increasing 
temperature (HORAI et al., 1970; MURASE and McBIRNEY, 1970). 

This effect will not introduce a substantial error in the 
estimate of the heat flux which we find to be 1.6x10 ' cal/cm sec 
using the "basalt" composition for the mantle; olivine would 
double this estimate. The former value is approximately 1/8 
that of the Pre-Cambrian regions of Earth (KAULA, 1968). For 

a Moon scaled Earth, the flux would be twice that actually 

measured for a "basalt" model and approximately equal for an 
ol i vi ne model . 

We wish to thank our colleaques C.W. Snyder, M. Neugebauer, 
J.H. Wolfe, and the Vela group at Los Alamos Scientific Labora- 
tory for supplying plasma data for testing the field line con- 
finement problem. Constructive comments by W.M. Kaula have 
substantially improved this paper. 
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APPENDIX 


The spherical harmonic expansions for the potential and 
magnetic field of the transverse electric part of the solar 
wind excitation are 

co 

0 = M-vHq Y rj t (-^—) (cose) 

- 6=1 

dp|( cos 9 ) 

w — 

p|( cos9 ) 

”'™s7FPn~ 


H r = H 0 sincp£ t(*+l ( c ° sS >- 

- 6=1 

The TE part of the magnetic field is related to the TE 
potential by 

H r = " (“S + k2 ) (rQ) » H 0 = TF ( r If) » 


K cp " pY?sTrfS' TF ( r If) 
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Thus the TE magnetic field in the lunar interior is 

dp|( cos9 ) 

W~ 

p|( cos 0 ) 
s i n 0 ' 

H r = Hq sincp X t(t+l) (fi) P{ (COS8) . 

1 = 1 
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FIGURE CAPTIONS 


Fig, 1 Power spectral density determinations for simultaneous 

magnetic field observations on the lunar surface (Apollo 
12) and in the solar wind near the Moon (Explorer 35). 
The y component is east and the z component north 
at the Apollo 12 site. Significant amplification is 
observed, increasing with frequency. Length of the 
time series is 2 hours. 

Fig. 2 Amplitude of the transfer function (ratio of lunar 
surface to free stream magnetic measurements) vs. 
frequency. The error bars for frequency are the 
windows defined by the lags in the autocorrelation 
calculation. The error bars in amplification are the 
one standard deviation limits determined from the means 
of 14 data spectra. Strong amplification is observed 
in the tangential (north and east) components of the 
field, while for the normal component the amplification 
follows the expected value of unity. 

Fig, 3 Distribution of the values of the differences of 
7\(f) from the mean °K at each frequency for all 
frequencies up to 0,035 Hz in units of the normalized 
standard deviation. 
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Fig. 4 The rms tangential lunar transfer function* A = 

[0.5(A^ + A^)] 2 and the normal transfer function A 
y z x 

as a function of frequency. The error bars, frequency 
windows, and directions are as shown in Fig. 2, The 
solid line is the value of ampl i f i cat ion calculated 
from the conductivity profile whose correspondi ng 
amplification is fitted to experimental values at 
frequencies of 0.83, 1.7, 5, 12, 17, 22, 25, and 35 
mill ihz. 

Fig. 5 Lunar bulk electrical conductivity profiles determined 
from the individual transfer functions A^ , A 2 , and 
7\ . The prominent rise of a from the surface inwards 
to r = 1500 km is apparent for all three cases as 
well as the subsequent decrease inwards to r = 1400 km 
followed by a more gradual rise. A tentative version 
of a lunar thermal profile is shown as the grey over- 
lay with temperatures indicated on the right hand margin. 
This profile is a fit of conductivities to a Nagata 
basalt in the mantle, an England olivine in the core, 
and the known subsurface temperature of -30°C. The 
insert is shown to suggest extreme values of the o's 
using the one standard deviation limits of the A's 
to calculate conductivities. 
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Fig, 6 Profiles of magnetic field amplitude and phase for the 
conductivity profile labelled in Fig, 5. The values 
are normalized to the transfer function at the lunar 
surface. The region of high conductivity centered 
at r = 1500 km (see Fig, 5) is responsible for strong 
attenuation and large phase shift at the higher 
frequencies. There is a correspondence between the 
selective attenuation shown here and the amplification 
vs, frequency curve of Fig. 4, as discussed in the 
text. Negative phase angle in this figure represents 
a time delay with respect to the driving function. 

Fig, 7 The experimental r.m.s. tangential lunar transfer 

function, S' , as a function of frequency. The error 
bars are the same as those in Figs. 2 and 4. The 
curve labelled "best fit" are the values of the 
amplification for the conductivity profile labelled 
7\ in Fig. 5. The "2 layer" curve was calculated for 
a model with a core of radius 1560 km and a constant 
conductivity of 7.6xl0 -if mhos/m, and a shell of zero 
conductivity. This value of the core conductivity 
provides a best fit to 7S for the given core radius. 

The DYAL and PARKIN (1971) amplification curve was 
calculated for a model with an inner core of radius 
1044 km and a conductivity of 1 0 ' mhos/meter, a 
middle layer extending to 1653 km with a conductivity 
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-4 . 

of 1.7x10 mhos/m and an outer layer of zero conductivity. 

Other amplification curves are based on conductivity 

models proposed by NESS (1969) and SILL (1971). Note 

that the "2 layer" curve does not duplicate the "S" 

nature of the experimental data. 

g. 8 Model Moon based upon the magnetometer data, the 
derived conductivity profiles, and the consistent 
fit of ct to rocky material. The anorthositic crust 
is taken from WOOD et al . (1970); it is not seen in 
the present analysis. 
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